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Abstract. We study the time behavior of entanglement between two noninteracting 
qubits each immersed in its own environment for two different non-Markovian 
conditions: a high-Q cavity slightly off-resonant with the qubit transition frequency and 
a nonperfect photonic band-gap, respectively. We find that revivals and retardation of 
entanglement loss may occur by adjusting the cavity-qubit detuning, in the first case, 
while partial entanglement trapping occurs in non-ideal photonic-band gap. 
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1. Introduction 



Real quantum systems unavoidably interact with their surroundings undergoing a 
consequent decoherence and entanglement loss [lj. It is known that two entangled 
qubits embedded in Markovian (memoryless) environments may become completely 
disentangled at a finite time, in spite of an exponential decay of the single qubit 
coherence [21 [3]. This entanglement sudden death (ESD, or early-stage disentanglement) 
which has been experimentally revealed [HE], puts a serious limit to the storage times of 
entanglement for practical purposes, e.g., for the realization of quantum memory banks 
[6]. A realistic quantum computer will probably have to take into account this quantum 
dynamical drawback. It is therefore important to study the possibile physical conditions 
where entanglement can be maintained. 

Entanglement losses during the evolution crucially depend on the particular noise 
acting on the system. Under Markovian-noise conditions the quantum process is 
typically irreversible. Differently, structured environments or strong coupling can 
give rise to non-Markovian noise (environment with memory) whose effects on the 
entanglement dynamics is currently subjected to investigation [6J. In this context, 
the cases of two noninteracting qubits embedded either in separated high-Q cavities 
[3 E] or in a common cavity [9] supporting a mode resonant with the qubit transition 
frequency have been analyzed, finding that revivals of the initial two-qubit entanglement 
can occur. When the qubits share a common environment, it has been also shown 
that entanglement can be preserved by means of quantum Zeno effect [10J. Moreover, 
entanglement trapping is achievable when two independent qubits are embedded in an 
ideal photonic-band gap (PBG) material (photonic crystal) [TT| [12]. 

The aim of this paper is therefore to deepen the analysis of entanglement evolution 
in non-Markovian environments, considering in particular two different effective spectral 
conditions of the environment-qubit system simulating, respectively, a cavity with a 
mode slightly off-resonant with the qubit transition frequency and a nonperfect photonic 
band-gap at the qubit transition frequency. In this paper we highlight the differences 
with previous studies and we discuss the optimal physical parameters for observing 
entanglement revivals and for preservation of entanglement. 



2. Model 

We consider a system composed by two parts S = A, B, each one consisting of a two-level 
system (qubit) S = A, B interacting with a reservoir R s = Ra, Rb- 

The single part S "qubit S +reservoir Rs" is described by the Hamiltonian 

H s = fiu a + d-- + ^2 ft [Ukblh + (gkfr+h + 9k°-bk) > (!) 
k 

where uq is the transition frequency and a± are the qubit raising and lowering operators, 
b\, b k are the creation and annihilation operators and g% the coupling constant of the 



Two-qubit entanglement dynamics for two different non-Markovian environments 3 



mode k with frequency lo^. When the environment is at zero temperature the single- 
qubit reduced density matrix p s {t) can be written, in the basis {|1), |0)}, as [1] 

/ pIMW)\ 2 p? (%(0 \ 



P b (t) 



(2) 



\ P&(0)<Z*(t) P^ (0) + pf 1 (0)(l-|g(t)| 2 ) J 

From the equation above it is readily seen that the single-qubit dynamics depends 
only on the function q(t) that in turn is determined by the reservoir spectral density. 
Indeed, q(t) obeys the differential equation q(t) = — f Q dtif(t — ti)g(ti), and the 
correlation function fit — t\) is related to the spectral density J(u) of the reservoir 
by fit — ti) = J duiJ{uj) exp[i(uo — uS)(t — ti)]. The solution of the associated algebraic 
equation for q(t), obtained trough Laplace transforms, is q(s) = q(0)/[s + f(s)], where 
q(s) and f(s) are the Laplace transforms of q(t) and fit — ti). 

A crucial quantity for our study is the explicit analytic expression of the two- 
qubit reduced density matrix at the time t. For the system of noninteracting qubits 
in separated environments here considered this can be obtained by a procedure 
based on the knowledge of the single-qubit dynamics [7J. In fact, given the time- 
dependent single-qubit density matrix elements as pf^it) = Ylw ^'(*)Pzz"(0)> Pfj'ft) = 
Emm'^jT'W' 9 ™^)) the time-dependent two-qubit density matrix elements are [7] 



P*JJ>(t)= E ^(*)^'(*)Phw(0), (3) 

11' mm' 



where i,j,l,m = 0,1. The two-qubit density matrix p(t) is thus obtained by means 
of Eq. ([3]) for an arbitrary two-qubit initial condition, its elements depending only on 
their initial values and on the function q(t) [3 [8]. In the following these density matrix 
elements will be meant in the standard computational basis E = {|1) = |11),|2) = 
|10),|3) = |01),|4) = |00)}. 



2.1. Initial states 

With regards to the initial state, we will limit our analysis to the case of initial pure 
Bell-like states 

|$) = a |oi) + /3e i<5 |10>, = a|00) + /3e i<5 |ll>, (4) 

with a, (3 real and a 2 + (3 2 = 1. For a = ±/3 = l/y/2 these states coincide with the Bell 
states. Bell-like states have the property that their resulting density matrix has an X 
structure (only diagonal and antidiagonal density matrix elements different from zero). 
Under our dynamical conditions, the X structure is maintained during the two-qubit 
evolution, so that the two-qubit density matrix at time t will also have a X structure. 



2.2. Concurrence 

In order to quantify the entanglement during the evolution of the bipartite system, we 
use the concurrence C [13J. The concurrence at the time t for an initial general X state 
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can be easily computed by exploiting the fact that the X structure is preserved here and 
by using Eqs. (j2J) and Q. The expression of the concurrence is given by [H 

C?(t) = 2max{0,K 1 (t),K 2 (t)}, 
Ki(t) = |g(t)| 2 {p 23 (0) - ^Pii(0)[p 4 4(0) + pii(0) 

x (1 - k(t)| 2 ) 2 + (p 22 (0) + P 33(0))(l - |g(t)| 2 )]- 
K 2 (t) = |g(t)| 2 [p 14 (0) - vW0) + pn(0)(l-|g(t)| 2 ) 

x v / P3 3 (0)+Pn(0)(l-|g(t)| 2 )] . (5) 

These formulas are quite general since their form does not explicitly depend on the 
particular choice of the environment, but only on the Hamiltonian model of Eq. (PQ) 
and on the chosen initial state. The explicit time dependence of concurrence depends 
on the explicit form of the function q(t) and thus it contains the information about 
the environment structure. In the following we shall consider two different environment 
structures with given spectral densities which in turn shall determine the explicit form 
ofg(t). 



3. Spectral density effect on the entanglement dynamics 

We shall now analyze the evolution of concurrence for two different spectral densities: a 
single Lorentzian simulating a cavity with a mode nonresonant with the qubit transition 
frequency and a nonperfect photonic band-gap at the qubit transition frequency. 



3.1. Off-resonant high-Q cavity 



As a first example we take the spectral distribution J(uj) of the electromagnetic field 
inside a high-Q cavity supporting a mode detuned by A from the qubit (two-level atom) 
transition frequency Uq, resulting from the combination of the environment spectrum 
and the system-environment coupling. It has the Lorentzian form [1] 

1 TA 2 

JM = 2^ -A-c) 2 + A 2 ' (6) 
where T is the qubit free space linewidth and A the spectral width of the coupling. 
The parameter A is then connected to the reservoir correlation time tb by the relation 
tb ~ A -1 . The relaxation timescale r# over which the state of the system changes 
is related to r by tr w T^ 1 . The correlation function corresponding to this J(uj) is 
f(t — £i) = ^exp[— (A — iA)(t — ti)}. Using this correlation function, the Laplace 



transform of q(t) is q(s) 
gives 

q(t) = e"- 



V 



s + 



I'A 



1 

2 s-(A-iA) 



and inverse Laplace transform finally 



cosh I — 



+ 



A -iA 



d 



sinh ( — 
2 



(7) 



where d = y/(X - iA) 2 - 2T\. 
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Figure 1. Nonresonant cavity: A = 0.1T. Concurrence as a function of the 
dimcnsionless quantities Tt starting from the initial EWL state /3*(0) = with 
a = l/\/3 for different values of detuning A: A = (solid curve), A = 2A (dotted 
curve), A = 5 A (long-short-dashed curve), A = 8A (long-short-short-dashed curve). 

In the resonant limit, A = 0, the correlation function has an exponential form with 
A the decay rate. In the analysis of the function q(t) of Eq. (J7|) with A = a weak 
(T < A/2) and a strong (r > A/2) coupling regime can be distinguished [H [15]. In 
the weak coupling regime the relaxation time is greater than the reservoir correlation 
time {tr > 2tb) and the behavior of q(t) is essentially a Markovian exponential decay 
controlled by T. In the strong coupling regime, the reservoir correlation time is greater 
or of the same order of the relaxation time (t>j < 2t#) and non-Markovian effects 
become relevant. Within this regime, q(t) presents oscillations describing a quasi- 
coherent exchange of energy between the qubit and the reservoir. 

3.1.1. Entanglement dynamics. We now investigate the entanglement dynamics of 
the two-qubit system considering in detail the effect of the detuning A, by using the 
concurrence obtained by Eq. ([5]) with q(t) given by Eq. (CO). The evolution of concurrence 
for various values of the detuning is shown in Fig. [1] when the cavity bandwidth A is 
smaller than the free-space atomic linewidth T (A = 0.1T) and the initial state is not 
maximally entangled, in particular p*(0) = | X I / )(^| with a = l/v^- In the resonant 
case it is found, as known, that the state suffers ESD [7j. However, from the plot one 
observes that increasing the detuning, the entanglement decay slows down. In particular, 
for A = 2A we have also revivals of entanglement after a finite period of time when the 
two qubits are not entangled. The phenomena of the slowing down of entanglement 
decay and of the entanglement revivals are a clear manifestation of the environment 
memory effects. 
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3.2. Nonperfect photonic band-gap 

As a second example of non-Markovian environment, we consider a spectral density of 
the form [16] 

T(u) _ J_ ( M ^x^_\ 

J[U) 2vr V>-^ ) 2 + A? {u-oo Q f + \l)' W 

which represents a Lorentzian with a dip used as a model to simulate the spontaneous 
decay of a qubit in a nonperfect photonic band gap (PBG). In Eq. (|Sj), Ai represents 
the bandwidth of the flat background continuum, A 2 the width of the gap, and T 1 and 
T 2 the strength of the background and the gap, respectively. The spectral density must 
be positive, this implies TiXf > I^Ag (condition for J{oj) to be positive at large to) and 
r\ > r 2 (condition for J (a;) to be positive at the center of resonance). Combining these 
two relations the condition I^Ai > r 2 A 2 must be satisfied, this being also the condition 
for a localized dip [IB] . In the case Y\ = T 2 the spectral density reduces exactly to zero at 
the center of the gap (u = Uq), leading to population trapping [IB] . For this form of J(uj), 
one obtains for the correlation function f(t — ti) = (r 1 Aie _Al ^' _ ' 1 ^ — r 2 A 2 e _A2 ( <_<1 ^) /2 
and the Laplace transform of q(t) becomes 

n( 8 ) = (A! + g )(A 2 + g ) 

qy ' s 3 + s 2 (Ax + A 2 ) + s(AxA 2 + A) + AiA 2 lY 1 ; 

where A = (FiAi — T 2 A 2 )/2 and = (I^ — T 2 )/2. By inverting the Laplace transform 
one finally obtains 

q(t) = T ^ + ^ (Al +/ 2) + A /V f , (10) 
^ {Ui — Uj){Ui — U k ) 

where i,j, k = 1, 2, 3 are all different indexes and Ui are the three solutions of the third 
degree equation appearing in the denominator of Eq. (jUJ) . The explicit expressions of 
are rather complex and we shall not report them here. 



3.2.1. Entanglement dynamics. We are now ready to analyze the entanglement 
dynamics of the two-qubit system in a nonperfect PBG as described by the spectral 
density of Eq. (jSJ), by using the concurrence of Eq. (j5J) with q(t) given by Eq. (TTOT) . In 
Fig. [2] we investigate the evolution of concurrence for Ai = 10A 2 = 50r! and varying 
T 2 . As initial state we consider the Bell state p*(0) = |$)($| (a = (3 = l/y/2). For 
T 2 = Ti, the spectral density goes to zero at the center of the gap and as a consequence 
we obtain entanglement trapping (similarly to what happens in Ref. [H]). The other 
curves are obtained for decreasing values of T 2 : T 2 = 2r 2 /3, T 2 = Ti/3 and T 2 = 0. 
Smaller values of T 2 correspond to a smaller dip of the spectral density at ui = ujq. In 
particular, for T 2 = 0, the shape of the spectral density is again a simple Lorentzian. In 
Fig. [2] the value of Ai is chosen so that we are in weak coupling regime (see Sec J3. II) and 
a Markovian decay occurs. Finally, the plot evidences how decreasing T 2 the trapping 
of entanglement is lost and the entanglement decay speeds up always more. 
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Figure 2. Nonperfect PBG case: Ai = IOA2 = 50Fi. Concurrence as a function 
of the dimensionless quantity Tit starting from the initial state p*(0) = |$)($| with 
a = (3 = l/y/2 for different values of P 2 : r2 = Ti (solid curve), r 2 = r 2 /3 (dotted 
curve), T 2 = 2Ti/3 (long-short-dashed curve), T 2 = (long-short-short-dashed curve). 



4. Conclusions 

In this paper we have extended previous analysis on entanglement dynamics of two 
noninteracting qubits embedded in bosonic environments at zero-temperature. We have 
examined two different spectral densities corresponding to two different environments: 
the first case considered is a Lorentzian spectrum representing a high-Q cavity out 
of resonance with the qubit transition frequency; the second, a nonperfect photonic- 
band gap. The first case has allowed to analyze the role of the cavity-qubit detuning 
comparing it to the known resonant case. In particular, as expected, an increase of the 
entanglement lifetime is observed when the detuning is increased. On the other hand, 
the second spectral density has permitted the study of entanglement dynamics when 
ideal conditions of the photonic-band gap (spectral density exactly equal to zero for a 
frequency equal to the qubit transition frequency) are not satisfied. In this case, it has 
been found that increasing the value of the spectral density in the central frequency 
(qubit transition frequency) entanglement trapping, which is expected for the ideal 
case, vanishes while entanglement decay speeds up. This study has thus enlarged 
the knowledge of entanglement evolution under different non-Markovian conditions, 
providing more hints for future investigations on this topic. 
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